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In  recent  years  the  Earth's  gravitational  field  has  been  determined  with 
continually  improving  accuracy,  using  hundreds  of  thousands  of  observations  of 
Earth  satellites,  chiefly  optical,  laser  and  Doppler,  together  with  surface 
gravimetry  and,  most  recently,  altimeter  measurements  from  the  Geos  3 satellite. 
The  geopotential  is  usually  expressed  as  a double  series  of  tesseral  harmonics, 
and  several  hundred  of  the  harmonic  coefficients  are  evaluated. 

‘Progress  in  this  work  during  the  1970s  is  briefly  outlined,  and  some  attempt 
is  made  to  assess  the  accuracy  of  current  geoid  maps  and  sets  of  harmonic 
coefficients,  as  exemplified  in  the  latest  models  derived  at  the  Goddard  Space 
Flight  Center.  The  harmonic  coefficients  of  order  14,  15  and  30  in  the  Goddard 
Earth  Model  I0B  are  compared  with  values  obtained  independently  by  analysis  of 
resonant  orbits:  the  results  suggest  that  the  values  in  GEM  10B  are  realistic 
for  these  orders,  and  presumably  others.  It  appears  that  the  accuracy  of  the 
geoid  maps  is  now  approaching  1 m . ^ ^ 

This  Report  is  based  on  a paper  of  the  same  title  presented  at  a Royal 
Society  Discussion  Meeting  on  "Satellite  Doppler  tracking  and  its  geodetic 
applications",  held  on  10  and  II  October  1978. 
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1 PREAMBLE 

A comprehensive  description  of  the  Earth's  gravity  field  cannot  be  fitted 
into  a short  paper.  Here  the  aim  is  less  ambitious  - to  give  some  idea  of  recent 
progress  in  the  evaluation  of  the  geopotential  and  to  offer  a few  fragmentary 
answers  to  the  question  "How  accurate  are  the  values  of  the  harmonic  coefficients 
in  recent  geopotential  models?" 

In  this  paper  'gravitational  field'  and  'geopotential'  refer  to  the  gravi- 
tational attraction  due  to  the  mass  of  the  Earth  and  atmosphere.  'Gravity  field' 
refers  to  gravity  as  measured  at  the  Earth's  surface,  including  the  effect  of  the 
Earth's  rotation. 

2 EXPRESSION  FOR  THE  GEOPOTENTIAL 

It  is  usual  to  express  the  Earth's  gravitational  potential  U at  an 
exterior  point  (r,  0,  X)  as  an  infinite  series  of  tesseral  harmonics  in  the  form* 


■■fit  ®‘-: 


(cos  0){C£m  cos  + sin  mX[N 


&m 


(1) 


1=2  m=0 


where  r is  the  distance  from  the  Earth's  centre,  9 is  co-latitude,  X is 

longitude  (positive  to  the  east),  GM  is  the  gravitational  constant  for  the 
3 2 

Earth  (398600  km  /s  ) , R is  the  Earth's  equatorial  radius  (6378.1  km), 

P™(cos  0)  is  the  associated  Legendre  function  of  order  ro  and  degree  £.  , and 
and  are  the  normalized  tesseral  harmonic  coefficients,  which  require 

to  be  evaluated.  For  m 1 the  normalizing  factor  N^  is  given  by 

n2  = 2(211  + 1)(£  - m)!  (r 

i,m  (i  + m)  ! 

2 

For  m = 0 , however,  N = 2£  + 1 . 

X/U 

Other  representations  of  the  geopotential  are  possible,  but  this  format 
has  proved  to  be  most  convenient  in  studies  using  satellites.  In  practice  the 
series  is  arbitrarily  truncated  at  the  order  and  degree  beyond  which  it  is  not 
thought  possible  to  determine  meaningful  values  of  the  harmonic  coefficients 

and  . In  recent  years  the  truncation  has  been  made  at  orders  between 

16  and  36,  but  larger  arrays  of  coefficients  may  be  used  in  future. 

Equation  (1)  gives  the  exterior  gravitational  attraction  due  to  the  mass 
of  the  Earth  and  atmosphere.  If  the  gravity  felt  at  the  Earth's  surface  is  to 
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be  calculated,  a 'centrifugal  potential'  Jr  w sin  0 has  to  be  added  to  U , 
where  co  is  the  Earth's  angular  velocity  (72.92115  * 10  ^ rad/s),  and  U should 
be  the  potential  of  the  Earth  excluding  the  atmosphere. 

To  obtain  a pictorial  expression  of  the  tesseral  harmonics  in  equation  (1), 
it  is  useful  to  think  of  the  suffix  m as  specifying  variations  from  one 
meridian  to  another.  A harmonic  of  order  m (for  any  degree  l >*  m)  exhibits 
m undulations  as  the  longitude  A increases  by  360°,  for  a fixed  latitude,  as 
shown  in  Fig  1 for  m = 15.  For  any  fixed  value  of  m , the  suffix  l determines 
the  variation  from  one  latitude  to  another.  A harmonic  of  degree  l and  order 
m has  (l  - m)  zeros  in  going  from  pole  to  pole  along  a fixed  longitude,  exclud- 
ing zeros  at  the  poles.  If  m = 0 , there  are  no  zeros  at  the  poles,  and  a 
section  of  the  Earth  through  the  poles  would  exhibit  fe  humps  in  360°  due  to 
the  (H,  0)  harmonic.  If  m(>0)  is  even,  there  are  zeros  at  the  poles  which  are 
maxima  or  minima,  and,  including  these  in  the  count,  a slice  through  the  poles 
would  show  (£  - m + 2)  humps  in  360°.  If  m is  odd,  the  polar  zeros  are  neither 
maxima  nor  minima  (being  points  of  inflection) , so  the  number  of  humps  is  reduced 
to  (1  - m + 1 ) . 

3 THE  SMITHSONIAN  STANDARD  EARTH  II 

Coefficients  of  some  tesseral  harmonics  in  the  Earth's  gravitational  field 

were  successfully  evaluated  in  the  1960s,  particularly  those  of  low  order,  but 

the  first  satisfactory  comprehensive  model  was  the  Smithsonian  Standard  Earth  II 
2 

published  in  1970  by  the  Smithsonian  Aetrophysical  Observatory,  Cambridge, 
Massachusetts.  The  Smithsonian  Standard  Earth  II,  which  was  a great  improve- 
ment on  its  predecessors,  relied  largely  on  100000  optical  observations  of 
satellites  from  Baker-Nunn  cameras,  with  an  observational  accuracy  of  about  10  m. 
The  expansion  of  the  geopotential  was  truncated  at  degree  and  order  16,  so  that 
there  were  about  250  geopotential  coefficients  to  evaluate.  The  orbital  pertur- 
bations caused  by  this  truncated  geopotential  were  calculated,  and  the  values  of 
the  geopotential  coefficients  and  station  coordinates  were  then  adjusted  so  that 
the  observations,  of  21  satellites  from  about  30  stations,  achieved  the  best 
possible  fit  to  these  perturbed  orbits,  and  also  satisfied  geometrical  constraints 
for  simultaneous  observations.  In  effect,  about  200000  equations  were  being  solved 
by  least  squares  for  more  than  300  unknowns,  namely  the  geopotential  coefficients 
and  the  station  coordinates. 

The  shape  of  the  gravity  field  given  by  the  Smithsonian  Standard  Earth  II 
is  shown  in  Fig  2.  This  map  gives  the  contours  of  the  geoid,  or  mean  sea-level 
surface,  relative  to  a reference  spheroid  having  a flattening  of  one  part  in 
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298.255,  and  the  heights  are  in  metres.  Stippled  areas  indicate  regions  where 
the  contour  heights  are  negative.  The  accuracy  is  of  order  5 m.  This  map  has 
the  advantage  of  being  easy  to  interpret:  it  indicates,  for  example,  that  if  you 
swam  along  the  equator  from  south  of  India,  where  there  is  a depression  of  about 
110  m,  to  north  of  New  Guinea,  where  there  is  a hump  of  about  80  m,  you  would  at 
the  end  of  your  marathon  swim  be  about  190  m further  from  the  Earth's  centre 
than  when  you  started,  though  you  never  go  uphill.  The  contours  define  the 
rather  peculiar  shape  taken  up  by  the  Earth's  sea-level  surface  in  response  to  the 
gravitational  pull  of  the  rather  peculiar  distribution  of  mass  inside  the  Earth. 
Apart  from  the  high  area  near  New  Guinea,  there  are  also  major  humps  near  Britain 
and  south  of  South  Africa,  both  about  60  m high.  The  main  depression  is  that 
south  of  India,  but  there  are  three  others  between  45  and  65  m deep,  south  of 
New  Zealand  and  near  Florida  and  California.  Since  the  pull  of  gravity  (if 
tides  and  transient  effects  are  averaged)  acts  in  a direction  perpendicular  to 
the  mean  sea- level  surface,  a geoid  map  like  that  of  Fig  2 gives  probably  the 
clearest  pictorial  representation  of  the  gravity  field.  In  terms  of  gravity 
anomalies,  that  is  the  difference  between  the  measured  acceleration  due  to 
gravity  and  a reference  field  symmetrical  about  the  equator,  the  picture  looks 
rather  different.  The  greatest  negative  anomaly,  about  -60  mgal,  is  south  of 
India,  and  the  largest  positive  anomaly,  about  40  mgal,  is  near  New  Guinea,  but 
the  other  main  features  are  more  numerous  and  differently  placed,  with  +33  mgal 
in  Alaska  and  +32  mgal  in  eastern  Europe,  for  example. 

4 THE  GEM  10  GEOID 

Since  1970,  many  new  gravitational  field  models  have  been  published,  includ- 
ing the  Smithsonian  Standard  Earths  III  and  IV  (Refs  3 & 4) , the  US  Department  of 
Defense  World  Geodetic  System  1972  (Ref  5),  the  European  models  GRIM  1 and  2 
(Refs  6 & 7),  and  the  series  of  Goddard  Earth  Models  developed  at  the  Goddard 
Space  Flight  Center,  which  have  appeared  in  pairs,  GEMs  1 and  2,  GEMs  3 and  4, 
and  so  on,  the  latest  being  GEMs  I0A  and  I 0B  (Ref  8).  The  observations  used  have 
gradually  been  extended  during  the  1970s.  As  well  as  the  photographic  observa- 
tions, large  numbers  of  Doppler  measurements  from  the  Navy  navigation  satellite 
system  have  been  brought  in,  together  with  laser  observations  of  ever-improving 
accuracy,  surface  gravimetry  of  steadily  increasing  coverage,  and  recently, 
radar  altimetry  from  Geos  3,  as  well  as  smaller  amounts  of  other  data.  The  geoid 
maps  derived  from  these  models  fortunately  all  look  fairly  similar,  and  GEM  10, 
shown  in  Fig  3,  is  probably  the  best  available  (in  the  absence  of  a GEM  10B  map). 
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The  GEM  10  model  uses  840000  observations,  including  213000  laser  ranges, 
150000  optical  observations  and  270000  US  Navy  Doppler  measurements.  The  most 
recent  laser  ranges  are  the  most  accurate  of  these.  The  solution  also  utilizes 
a set  of  1654  equal-area  5°  surface  gravity  measurements*^,  which,  although 
worldwide,  are  weak  in  the  southern  oceans.  In  all,  592  harmonic  coefficients 

m 

are  evaluated  in  GEM  10,  and  the  geopotential  is  complete  to  degree  and  order  22. 
Fig  3 shows  the  same  main  features  as  Fig  2,  although  the  world  has  been  split  at 
longitude  0°  rather  than  180°.  The  depression  south  of  India  is  now  105  m deep 
as  against  113  m in  the  earlier  model,  but  the  other  three  dips  are  quite  simi- 
lar, 59  m as  against  61  m south  of  New  Zealand,  46  m as  against  45  m off  California 
and  52  m as  against  50  m off  Florida.  The  New  Guinea  hump  is  73  m as  against  81  m. 
The  hump  near  Britain  is  61  m,  the  same,  and  that  south  of  Africa  is  48  m as 
against  56  m. 

GEM  10  is  probably  accurate  to  I or  2 m,  if  you  accept  that  it  inevitably 
irons  out  any  fine  detail,  because  even  the  22nd  harmonic  has  a semi-wavelength 
of  8°,  or  900  km.  The  fine  detail  can  best  be  appreciated  by  looking  at  the 
preliminary  maps  of  the  ocean  surface  obtained  from  altimeter  measurements  by  the 
Geos  3 satellite  (1975-27A):  one  of  these  (from  Ref  11)  is  shown  in  Fig  4.  Here 
the  detail  is  on  a scale  of  a few  kilometres,  and  the  relative  accuracy  over 
small  areas  should  be  excellent,  though  the  absolute  accuracy  over  large  areas 
must  be  treated  with  caution.  The  original  maps  are  at  contour  intervals  of  1 m, 
but  for  clarity  only  the  contours  at  5 m intervals  are  shown  in  Fig  4.  The  con- 
tours in  this  section  of  the  North  Atlantic  to  the  west  of  Europe  differ  in 
detail  from  those  of  GEM  10,  being  much  more  tortuous,  but  the  general  trends  are 
similar.  At  a latitude  of  60°  N there  is  an  increase  of  about  13  m between 
longitude  0 and  30°  W in  the  GEM  10  map  (Fig  3),  and  15  m in  Fig  4.  At  longitude 
30°  W there  is  a decrease  of  only  about  5 m in  geoid  height  between  latitudes 
60°  N and  40°  N in  GEM  10,  and  the  value  is  the  same  in  Fig  4.  The  geoid  in  this 
area  is  very  flat  compared  with  many  other  areas  of  the  world. 

GEM  10  did  not  use  altimetry  data,  so  Fig  4 is  a nearly  independent  test 
which  confirms  the  accuracy  of  the  GEM  10  geoid.  Current  and  future  GEM  models 
include  altimetry  data,  which  naturally  tend  to  dominate  the  solutions  for  har- 
monics of  order  higher  than  about  30. 

5 ACCURACY  OF  THE  HARMONIC  COEFFICIENTS 

5. 1 The  test  of  resonance 

Although  the  GEM  geoid  contours  may  be  accurate  to  about  1 m (excluding 
fine  detail),  the  accuracy  of  the  592  harmonic  coefficients  is  still  questionable. 
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With  more  than  a million  equations  to  be  solved  for  more  than  600  unknowns. 
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are  bound  to  be  some  very  high  correlations;  or,  to  put  it  in  another  way,  very 
different  sets  of  coefficients  might  lead  to  nearly  the  same  geoid.  We  need  a 


correct  set  of  coefficients  not  only  to  give  the  correct  gravity  field,  but  also 
because  they  would  provide  a strong  indication  of  the  mass  distribution  in  the 
Earth's  interior,  supplying  a criterion  for  judging  between  existing  theories  of 


the  lithosphere  and  upper  mantle  of  the  Earth. 


In  the  early  1970s  the  coefficients  of  order  higher  than  about  10  were  not 
very  reliable,  apart  from  a few  for  which  there  were  orbits  in  shallow  resonance; 
but  recent  models  show  great  improvement.  The  accuracy  of  some  high-order 
ceofficients  can  be  indeper lently  checked  by  analysis  of  orbits  that  pass  through 
resonance  with  the  Earth's  gravitational  field.  'Resonance'  occurs  when,  after  a 
certain  number  of  revolutions,  the  satellite  repeats  its  track  over  the  Earth. 

For  example,  if  the  orbital  period  is  such  that  the  Earth  spins  through  exactly 
24°  relative  to  the  orbital  plane  between  one  equator  crossing  and  the  next, 
successive  tracks  of  the  satellite  over  the  Earth  will  be  24°  further  west.  After 
15  orbits  the  ground  track  will  have  moved  360°  and  will  then  repeat  itself. 

This  is  15th-order  resonance,  and  when  it  occurs  the  perturbing  effects  of  15th- 
order  harmonics  in  the  geopotential  build  up  day  after  day  until  there  is  quite 
a large  change  in  some  orbital  parameters,  particularly  the  inclination  of  the 
orbit  to  the  equator.  Accurate  measurement  of  the  change  in  inclination  will 
give  the  value  of  a linear  sum  of  harmonics  of  15th  order  and  odd  degree,  a 
'lumped  harmonic',  as  it  is  usually  called.  These  resonances  occur  as  the 
orbits  slowly  contract  under  the  influence  of  air  drag,  and  if  the  orbit  is  con- 
tracting slowly  enough,  the  change  in  inclination  can  be  very  accurately  fitted 

12 

with  the  appropriate  theoretical  curve.  Fig  5 shows  the  resonant  variation 
of  the  satellite  1971-54A  between  November  1972  and  September  1974.  The  inclina- 
tion decreased  by  about  0.04°,  equivalent  to  5 km  on  the  Earth's  surface;  meas- 
uring this  major  orbital  change  offers  the  opportunity  of  accurately  determining 
lumped  15th-order  harmonics. 


By  analysing  a number  of  deep  resonant  orbits  of  this  type,  at  different 

. . . . . . 12  13 

inclinations,  values  of  individual  !5th-order  coefficients  were  obtained  ’ 

. 14 

in  1975;  more  recently  , values  of  individual  14th-order  coefficients  of  degree 
14  to  22  have  been  derived.  The  values  of  the  14th  and  15th-order  coefficients 
in  recent  geopotential  models  will  now  be  compared  with  these  independent  results. 

5.2  14th-order  coefficients 


Fig  6 shows  the  14th-order  S coefficients  of  degree  14  to  22,  as  obtained 
in  Ref  14:  the  values  are  shown  as  circles  with  error  bars  of  length  2 sd. 
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The  values  are  arranged  so  that  they  increase  from  left  to  right:  the  degree 
£ of  the  coefficients,  marked  at  the  top,  is  consequently  in  a rather  random 
order. 

Also  shown  in  Fig  6 are  the  values  from  three  recent  comprehensive  geo- 

4 . . 

potential  models.  The  x signs  are  values  from  the  Smithsonian  Standard  Earth 
IV. 3;  the  + signs  are  values^  from  the  European  model  GRIM  2;  the  diamonds  indi- 
cate values  from  the  latest  available  GEM  model.  This  model  is  not  GEM  10,  for 
which  the  geoid  map  was  shown  in  Fig  3,  but  GEM  10B  (Ref  8),  which  goes  up  to  degree 
and  order  36,  incorporates  the  results  of  700  passes  of  altimeter  data  from 
Geos  3,  and  is  considerably  better  than  GEM  10. 

For  most  values  of  £ there  is  good  agreement  between  the  values  of  S^  ^ 
in  Fig  6 from  resonance,  and  those  in  the  geopotential  models.  GEM  10B  agrees 
particularly  well  for  £ = 14,  15  and  16,  and  differs  seriously  only  for  £ = 21. 

Fig  7 shows  the  corresponding  comparison  for  Here  the  agreement 

with  GEM  10B  and  SSE  IV. 3 is  very  good  for  £ = 15  and  17,  but  the  even-degree 
coefficients  are  less  satisfactory.  Of  course,  there  is  no  certainty  that  the 
values  from  resonance  are  completely  reliable;  there  is,  however,  no  value  of  £ 
for  which  all  three  models  agree  on  a value  of  ^ different  from  the  value 
obtained  from  resonance  analysis,  so  it  is  very  probable  that  the  values  from 
resonance  are  reliable,  although  their  accuracy  obviously  needs  to  be  improved 
by  analysing  further  orbits. 

It  is  not  too  surprising  that  the  l4th-order  coefficients  in  the  compre- 
hensive solutions  are  fairly  accurate,  because  all  the  solutions  incorporate 
results  from  satellites  that  are  quite  close  to  14th-order  resonance,  and  there- 
fore suffer  relatively  large  perturbations  due  to  the  14th-order  terms.  A more 
severe  test  is  the  15th-order  coefficients,  which  we  now  examine. 

5 . 3 15th-order  coefficients 

12  13 

In  Fig  8 the  15th-order  S coefficients  derived  from  resonance  analysis  ’ 
are  shown  as  circles  with  error  bars  of  length  2 sd.  For  odd-degree  coefficients 
the  resonant  solution  was  carried  up  to  £ = 33,  and  all  the  values  are  shown; 
but  for  the  even-degree  coefficients,  only  £ = 16  and  18  are  given  in  Fig  8, 
because  the  two  higher-degree  coefficients  evaluated  (£  = 20  and  22)  were  less 
accurate  and  less  reliable.  Again  the  values  are  arranged  so  that  they  increase 
from  left  to  right.  The  values  from  GEM  10B  and  SSE  IV. 3 are  shown  for  compari- 
son. GRIM  2 is  excluded  because  it  uses  the  resonance  results. 
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Fig  8 shows  chat,  for  five  of  the  twelve  values  of  £ , the  values  of 
S . from  GEM  10B  agree  well  with  those  from  resonance  analysis.  The  agreement 

j 1 J 

is  best  for  £ = 15  and  17,  and  perhaps  surprisingly,  for  £ = 31  and  33. 
Certainly  the  latest  GEM  model  comorms  to  the  resonance  values  much  better  than 
the  GEM  models  which  were  current  in  1975,  when  the  resonance  analysis  was  pub- 
lished. In  GEM  10B  the  coefficients  of  order  higher  than  30  were  determined 
solely  from  the  altimeter  data.  Perhaps  this  is  why  the  values  of  ^ for 

£ = 31  and  33  agree  so  well  with  the  values  for  resonance. 

Fig  9 shows  a similar  diagram  for  C . These  results  are  better  than 

j 1 J 

those  in  Fig  8.  The  agreement  with  GEM  10B  is  good  for  £=  15,  16,  17,  19,  23, 

25,  29  and  33.  It  is  probably  fair  to  conclude  that  most  of  the  coefficients  in 

-9 

GEM  10B  up  to  order  15,  and  possibly  beyond,  are  accurate  to  within  ±5  x 10 
5 . 4 Higher-order  coefficients 

The  GEM  10B  model  includes  harmonics  to  order  and  degree  36:  is  it  possible 
to  test  the  accuracy  of  any  of  its  higher-order  coefficients  by  resonance 
analysis?  In  principle,  the  answer  is  'yes'»  It  should  be  possible  to  assess 
the  accuracy  of  29th  and  31st-order  harmonics  using  29:2  and  31:2  resonance.  The 
first  analysis  of  29:2  resonance  was  made  in  1976  by  Doreen  Walker,  who  succeeded 
in  obtaining  a lumped  29th-order  harmonic  coefficient  from  analysis  of  the  orbit 
of  Ariel  1 (1962ol)  at  29:2  resonance'3.  Unfortunately,  it  is  not  possible  to 
evaluate  for  comparison  the  corresponding  lumped  harmonic  from  GEM  10B,  because 
coefficients  up  to  degree  45  (or  more)  contribute  to  the  lumped  harmonic,  and 
GEM  10B  only  goes  to  degree  36. 

Analysis  of  the  31:2  resonance  is  more  difficult,  because  of  the  high  drag 
which  is  inevitable.  The  first  lumped  values  were  obtained  by  Hiller  and 
King-Hele*8  from  the  orbit  of  Proton  4 (1968-103A),  and  better  values  have 
recently  been  obtained'7  from  Skylab  1 rocket  (1973-27B).  The  equation  obtained 
for  the  S-coeff icients  from  Skylab  1 rocket  is: 

0.0.1*32,31  " °-08^34,31  + °-31°S6,31  " °-701?38,31 

+ ^40,31  " 0,799842,31  + °-084^44, 3J  + 0-508S46>31 

" °-372S48,3.  ~ °-192S50,31  + = (-  13.5  ± 2.7)  x 

Although  the  lumped  value  is  of  good  accuracy,  it  is  not  possible 
value  from  GEM  10B  for  comparison,  because  values  of  the  harmonic 


to  obtain  a 
coefficients 
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up  to  degree  50  or  beyond  would  be  required.  The  equation  for  the  lumped 
C-coeff icient  has  the  same  numerical  factors  on  the  left  hand  side,  and  the 

_9 

value  on  the  right-hand  side  is  (9.1  ± 4.2)  * 10 

Although  the  29th  and  31st-order  coefficients  in  GEM  10B  cannot  yet  be- 

tested  by  this  method,  the  30th-order  harmonics  are  open  to  verification.  In 

12. 

the  analysis  of  15th-order  resonance  it  was  possible  to  obtain  good  values  of 
lumped  30th-order  harmonics  from  four  satellites,  and  for  one  of  these,  the 
satellite  1971-54A  shown  in  Fig  5,  the  multiplying  coefficients  in  the  lumped 
harmonic  decrease  rapidly  as  the  degree  increases,  so  that  only  the  coefficients 
of  degree  30,  32,  34  and  36  are  needed.  The  lumped  S-coeff icient  from  1971 —54 A 
is  given  by: 


^30  = S30,30  + °'428S32,30  + °’211S34,30  + °'°97S36,30  + 0(0'03S38  ,30) 


= (15.3  ± 1.3)  x 10 
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For  the  C-coeff icient  the  corresponding  value  is: 


C^q2  = (-  10.3  ± 1.5)  x 10  9 . 

The  values  of  the  same  lumped  harmonic  coefficients  from  GEM  10B  are: 

-0  2 -9  -02  -9 

S3’  = 11.2  x 10  ’ and  CiJ*  = - 8.4  x 10 


So,  as  shown  in  Fig  10,  the  GEM  10B  values  are  in  the  same  direction  as  those 
from  resonance,  and  within  25%  of  the  numerical  value.  Unless  the  agreement  is 
just  luck,  it  seems  that  the  30th-order  coefficients  in  GEM  10B  may  be  quite 
realistic . 

5.5  Conclusions 

The  accuracy  of  the  geopotential  models  has  improved  greatly  in  recent 
years.  In  the  early  1970s  the  values  of  many  of  the  harmonics  of  order  greater 
than  10  were  rather  fictional,  and  the  15th-order  coefficients  determined  from 
resonance  differed  widely  from  those  in  the  models.  But  now  there  is  good  agree- 
ment with  GEM  10B  for  order  14  and  15,  and  it  is  probably  fair  to  conclude  that 

...  -9 

the  coefficients  m GEM  10B  should  be  accurate  to  ±5  x 10  (equivalent  to  about 

±20%  for  order  15),  for  degree  and  order  up  to  15,  and  perhaps  up  to  order  and 
degree  30  or  more,  if  the  good  agreement  of  the  one  result  available  for  compari- 
son is  not  fortuitous. 
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VALUES  OF  GM 


In  addition  to  the  values  of  the  harmonic  coefficients,  the  constant  GM 

3 2 

in  equation  (1)  needs  to  be  evaluated.  The  value  is  close  to  398600  km  / s , 

and  until  recently  the  best  method  of  measuring  this  constant  was  from  the 

trajectories  of  space  vehicles.  Values  of  (GM  - 398600)  obtained  from 

Mariners  9 and  10  and  Vikings  1 and  2 respectively  have  been  0.66  ± 0.06  (Ref  18), 

0.45  ±0.2  (Ref  19),  0.40  ± 0.2  and  0.60  +0.2  (Ref  20).  The  unweighted  average 

is  0.53.  Values  from  lunar  laser  ranging  are  0. 48  ± 0.1  (Ref  21)  and 

0. 52  ± 0.03  (Ref  22).  The  most  recent  and  probably  the  most  accurate  value,  from 

laser  tracking  of  near-Earth  satellites,  particularly  Lageos  (1976-39A),  is 

3 2 

GM  = 398600.44  + 0.02  km  /s  (Ref  23).  All  the  values  quoted  above  are  from 
Ref  23  where  the  previous  values  have  been  adjusted  to  a consistent  value  for 
the  speed  of  light,  namely  299792.458  km/s. 

The  value  of  the  Earth's  equatorial  radius,  R , from  GEM  10B  is  63781 39 ± 1 m. 
7 FUTURE  TRENDS 

During  the  1970s  the  accuracy  of  the  models  of  the  Earth's  gravity  field 
has  steadily  improved,  largely  as  a result  of  the  improving  accuracy  of  the  satel- 
lite laser  ranging  measurements,  now  supplemented  by  altimeter  data.  In  1970  the 
geoid  accuracy  was  5-10  m.  Today  the  figure  is  near  1 m,  and  the  values  of  indi- 
vidual harmonics  are  quite  realistic  up  to  order  15  and  possibly  order  30.  This 
continuing  advance  in  geodesy  contrasts  with  what  may  seem  to  be  comparative 
stagnation  in  solid-Earth  geophysics,  where  the  simplistic  concept  of  plate 
tectonics  provided  valuable  insight  in  the  late  1960s,  but  has  tended  to  harden 
into  dogma  without  maturing,  so  that,  for  example.  Earth  movements  are  often 
still  modelled  using  rigid  plates  of  constant  thickness. 

The  advances  in  geodesy  are  likely  to  continue.  There  are  many  new  methods 

24 

of  measurement  that  may  help,  for  example  satellite-to-satellite  Doppler  track- 
ing and  ranging.  But  progress  is  fairly  certain  without  any  new  techniques, 
merely  from  the  accumulation  of  mere  accurate  data  on  laser  ranging  to  satellites 
from  the  ground  and  a denser  mesh  of  more  accurate  altimeter  measurements.  In 
the  next  Goddard  Earth  model  it  is  probable  that  harmonics  up  to  order  and  degree 
180  will  be  evaluated,  so  that  there  will  be  about  30000  geopotential  coefficients 
to  determine  from  millions  of  observations.  In  this  180  x 180  field  the  coeffi- 
cients of  order  greater  than  about  30  will  presumably  be  determined  mainly  from 
altimeter  data.  A geopotential  field  complete  to  order  and  degree  180  may  seem 
rather  staggering,  but  a model  with  that  degree  of  detail  is  obviously  needed 
because  even  the  highest  harmonics,  of  order  180,  have  a wavelength  of  200  km. 
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and  the  altimeter  data  show  much  finer  detail.  So  the  future  promises  more 
accurate  and  more  detailed  geoid  maps,  and  a much  fuller  array  of  geopotential 
coefficients  of  improving  accuracy. 
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